Genome Evolution Research Articles

A high-quality chromosome-scale genome assembly of blood orange, an important pigmented sweet orange variety

Blood orange (BO) is a rare red-fleshed sweet orange (SWO) with a high anthocyanin content and is associated with numerous health-related benefits. Here, we reported a high-quality chromosome-scale genome assembly for Neixiu (NX) BO, reaching 336.63 Mb in length with contig and scaffold N50 values of 30.6 Mb. Furthermore, 96% of the assembled sequences were successfully anchored to 9 pseudo-chromosomes. The genome assembly also revealed the presence of 37.87% transposon elements and 7.64% tandem repeats, and the annotation of 30,395 protein-coding genes. A high level of genome synteny was observed between BO and SWO, further supporting their genetic similarity. The speciation event that gave rise to the Citrus species predated the duplication event found within them. The genome-wide variation between NX and SWO was also compared. This first high-quality BO genome will serve as a fundamental basis for future studies on functional genomics and genome evolution.

Teaching transposon classification as a means to crowd source the curation of repeat annotation – a tardigrade perspective

BackgroundThe advancement of sequencing technologies results in the rapid release of hundreds of new genome assemblies a year providing unprecedented resources for the study of genome evolution. Within this context, the significance of in-depth analyses of repetitive elements, transposable elements (TEs) in particular, is increasingly recognized in understanding genome evolution. Despite the plethora of available bioinformatic tools for identifying and annotating TEs, the phylogenetic distance of the target species from a curated and classified database of repetitive element sequences constrains any automated annotation effort. Moreover, manual curation of raw repeat libraries is deemed essential due to the frequent incompleteness of automatically generated consensus sequences.ResultsHere, we present an example of a crowd-sourcing effort aimed at curating and annotating TE libraries of two non-model species built around a collaborative, peer-reviewed teaching process. Manual curation and classification are time-consuming processes that offer limited short-term academic rewards and are typically confined to a few research groups where methods are taught through hands-on experience. Crowd-sourcing efforts could therefore offer a significant opportunity to bridge the gap between learning the methods of curation effectively and empowering the scientific community with high-quality, reusable repeat libraries.ConclusionsThe collaborative manual curation of TEs from two tardigrade species, for which there were no TE libraries available, resulted in the successful characterization of hundreds of new and diverse TEs in a reasonable time frame. Our crowd-sourcing setting can be used as a teaching reference guide for similar projects: A hidden treasure awaits discovery within non-model organisms.

Redefining the bacteriophage mv4 site-specific recombination system and the sequence specificity of its attB and core-attP sites.

Through their involvement in the integration and excision of a large number of mobile genetic elements, such as phages and integrative and conjugative elements (ICEs), site-specific recombination systems based on heterobivalent tyrosine recombinases play a major role in genome dynamics and evolution. However, despite hundreds of these systems having been identified in genome databases, very few have been described in detail, with none from phages that infect Bacillota (formerly Firmicutes). In this study, we reanalyzed the recombination module of Lactobacillus delbrueckii subsp. bulgaricus phage mv4, previously considered atypical compared with classical systems. Our results reveal that mv4 integrase is a 369 aa protein with all the structural hallmarks of recombinases from the Tn916 family and that it cooperatively interacts with its recombination sites. Using randomized DNA libraries, NGS sequencing, and other molecular approaches, we show that the 21-bp core-attP and attB sites have structural similarities to classical systems only if considering the nucleotide degeneracy, with two 7-bp inverted regions corresponding to mv4Int core-binding sites surrounding a 7-bp strand-exchange region. We also examined the different compositional constraints in the core-binding regions, which define the sequence space of permissible recombination sites.

Chromosome-level genome assembly of Pedicularis kansuensis illuminates genome evolution of facultative parasitic plant.

Parasitic plants have a heterotrophic lifestyle, in which they withdraw all or part of their nutrients from their host through the haustorium. Despite the release of many draft genomes of parasitic plants, the genome evolution related to the parasitism feature of facultative parasites remains largely unknown. In this study, we present a high-quality chromosomal-level genome assembly for the facultative parasite Pedicularis kansuensis (Orobanchaceae), which invades both legume and grass host species in degraded grasslands on the Qinghai-Tibet Plateau. This species has the largest genome size compared with other parasitic species, and expansions of long terminal repeat retrotransposons accounting for 62.37% of the assembly greatly contributed to the genome size expansion of this species. A total of 42,782 genes were annotated, and the patterns of gene loss in P. kansuensis differed from other parasitic species. We also found many mobile mRNAs between P. kansuensis and one of its host species, but these mobile mRNAs could not compensate for the functional losses of missing genes in P. kansuensis. In addition, we identified nine horizontal gene transfer (HGT) events from rosids and monocots, as well as one single-gene duplication events from HGT genes, which differ distinctly from that of other parasitic species. Furthermore, we found evidence for HGT through transferring genomic fragments from phylogenetically remote host species. Taken together, these findings provide genomic insights into the evolution of facultative parasites and broaden our understanding of the diversified genome evolution in parasitic plants and the molecular mechanisms of plant parasitism.

Update on feline calicivirus: viral evolution, pathogenesis, epidemiology, prevention and control

Feline calicivirus (FCV) is a prevalent and impactful viral pathogen affecting domestic cats. As an RNA virus, FCV exhibits high mutability and genetic plasticity, enabling its persistence within cat populations. Viral genetic diversity is associated with a broad spectrum of clinical manifestations, ranging from asymptomatic infections and mild oral and upper respiratory tract diseases to the potential development of virulent systemic, and even fatal conditions. This diversity poses distinctive challenges in diagnosis, treatment, and prevention of diseases caused by FCV. Over the past four decades, research has significantly deepened understanding of this pathogen, with an emphasis on molecular biology, evolutionary dynamics, vaccine development, and disease management strategies. This review discusses various facets of FCV, including its genomic structure, evolution, innate immunity, pathogenesis, epidemiology, and approaches to disease management. FCV remains a complex and evolving concern in feline health, requiring continuous research to enhance understanding of its genetic diversity, to improve vaccine efficacy, and to explore novel treatment options.

De novo transcriptome assembly and genome annotation of the fat-tailed dunnart (Sminthopsis crassicaudata)

Marsupials exhibit distinctive modes of reproduction and early development that set them apart from their eutherian counterparts and render them invaluable for comparative studies. However, marsupial genomic resources still lag far behind those of eutherian mammals. We present a series of novel genomic resources for the fat-tailed dunnart (Sminthopsis crassicaudata), a mouse-like marsupial that, due to its ease of husbandry and ex-utero development, is emerging as a laboratory model. We constructed a highly representative multi-tissue de novo transcriptome assembly of dunnart RNA-seq reads spanning 12 tissues. The transcriptome includes 2,093,982 assembled transcripts and has a mammalian transcriptome BUSCO completeness score of 93.3%, the highest amongst currently published marsupial transcriptomes. This global transcriptome, along with ab initio predictions, supported annotation of the existing dunnart genome, revealing 21,622 protein-coding genes. Altogether, these resources will enable wider use of the dunnart as a model marsupial and deepen our understanding of mammalian genome evolution.

Abstract PO5-16-01: Clinicogenomic analysis reveals genomic associations of brain metastatic tropism in breast cancer

Abstract Introduction: Recent progress in imaging and therapeutics has led to improvements in the treatment of advanced breast cancer. However, the incidence of brain metastasis (BM) in breast cancer is increasing and continues to be associated with poor prognosis and significant morbidity. Herein, we describe genomic enrichments associated with brain tropism from a large institutional clinico-genomic cohort. Methods: Patients who underwent tumor and matched normal sequencing utilizing hybridization-capture based targeted exonic sequencing (MSK-IMPACT) from April 2014 to March 2023 were included in this analysis; this cohort included 6,377 samples from 5,135 patients. We analyzed genomic data to inform the full spectrum of somatic and germline mutations, copy number changes and structural variants enriched in samples in resected BM specimens (n = 179), compared to primary tumors (n = 2980), stratified by receptor status. Enrichment analysis was conducted using a Fisher exact test for each respective receptor status. Similar analyses were conducted between extracranial metastatic (ECM) sites (n = 3186) and primary tumors, to differentiate the unique genomic features driving the emergence of BM as compared to general metastatic propensity. Multiple hypothesis testing was performed using Benjamini-Hochberg. Further genomic analysis, including allele specific copy number changes and mutational signatures will be reported at the 2023 SABCS meeting. Results: Oncogenic alterations in several genes were found to be enriched in BM samples compared to primary tumor, with an odds ratio exceeding that of the corresponding ECM comparison. For example, in the HR+/HER2- subset (4329 samples from 3476 patients), pathogenic variants in the following genes were enriched in BM:Primary analysis: NF1 (OR 6.4 [2.9 – 13.7], q = 0.001, 15.5% of BM), MYC (OR 3.9 [2.0 – 7.4], q = 0.001, 24.4% of BM), AGO2 (5.2 [2.4 – 10.9], q= 0.002, 15.5% of BM), CDKN2A (OR 3.9 [1.6 – 9.4], q = 0.06, 10.3% of BM), ERBB2 (OR 3.6 [1.5 – 8.6], q = 0.06, 10.3% of BM), RB1 (OR 4.2 [1.5 – 11.4], q = 0.06, 8.6% of BM), TP53 (OR 2.0 [1.1 – 3.5], q = 0.07, 40% of BM). In contrast, ESR1 variants were enriched in BM:Primary (OR 14.0 [6.7 – 28.2], but at the same level as ECM:Primary (OR 12.3 [8.7 – 17.7]). When considering HER2+ tumors (919 samples from 739 patients), only NKX.3 was enriched in the BM:Primary comparison (OR 11.5 [3.2 – 42.3], q = 0.004, 11.6% of BM). Although no TNBC samples met significance thresholds after adjustment for multiple hypothesis testing, CDH1, NTRK1, EGFR, and TP53 emerged as putatively significant. Conclusions: In a large cohort of genomically-profiled breast cancer samples, we find recurrent involvement of genes involved in cell cycle regulation (RB1, CDKN2A, TP53); these have been found to be enriched in brain metastases from other cancer types and are implicated in resistance to common lineage-directed therapies. Our approach also uncovered several targets for which therapies are actively in development (ERBB2, MYC). Further study of genomic evolution of breast cancer in the context of brain metastasis and therapeutic resistance will uncover conserved site-specific genomic changes and will identify novel, rational therapeutic targets. Citation Format: Anton Safonov, Deborah Smith, Emanuela Ferraro, Junchao Shen, Ishaani Khatri, Rahul Kumar, Julia An, Mark Robson, Sarat Chandarlapaty, Adrienne Boire, Nikolaus Schultz, Nelson Moss, Luke Pike, Pedram Razavi. Clinicogenomic analysis reveals genomic associations of brain metastatic tropism in breast cancer [abstract]. In: Proceedings of the 2023 San Antonio Breast Cancer Symposium; 2023 Dec 5-9; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2024;84(9 Suppl):Abstract nr PO5-16-01.

Abstract PO2-04-06: Investigation of the genomic evolution of HER2-positive breast cancer following progression on dual HER2-targetted therapy with Trastuzumab and Tucatinib

Abstract Background: Advanced HER2-positive breast cancer remains an incurable and highly morbid condition. One of the key challenges has been treatment of central nervous system (CNS) metastases, commonly seen in those with HER2-positive disease. Tucatinib is a recently approved treatment HER2-targetting tyrosine kinase inhibitor with proven CNS activity. There has been little analysis of potential genomic mechanisms of drug resistance. Method: We retrospectively identified patients that were treated with tucatinib, capecitabine and trastuzumab. We collated baseline clinical data and treatment history and all tumour and plasma samples available were obtained for testing. DNA was extracted from plasma samples using the QIAamp DNA Mini Kit (Qiagen) on 4ml plasma samples. Targeted next generation sequencing (NGS) of serial tumour and plasma samples was performed using a custom hybrid capture assay that covered 180 genes known to be recurrently mutated in breast cancer focusing on the known genomic landscape of metastatic disease as well as the inclusion of specific actionable targets. Survival analyses was performed using Cox regression models. Results: We identified 11 patients at our institution. Plasma samples were available for all 11 patients with a total of 107 serial samples across multiple time points, with a median of 8 (range 2-19) samples per patient. All 11 patients had baseline and End of Treatment (EOT) samples. In addition, 6/11 (54%) patients had tumour sequenced on the same NGS platform prior to commencement of therapy, with 12 individual tumour samples tested. The median PFS of our cohort was 9.6 months (4.0 - NR) with a median OS of 36.2 months (10.4-NR). 7/11 (64%) patients had CNS disease at baseline with their median PFS being 13.5 months (3.9-24.4) and median OS of 36.2months (10.4-NR). Overall, we observed an average of 18.4 (4-35) mutations per patient in tumour, 5.1 (1-15) in baseline plasma and 5.8 (2-14) in EOT plasma. The most frequent mutations seen in baseline tumour samples were ERBB2 (92%), KMT2C (67%) and NCOR1 (58%) compared with KMT2C (64%), NOTCH4 (36%) and PIK3CA (27%) in baseline plasma. Patients with a lower number of ctDNA mutations at baseline showed a trend towards improved survival (PFS 14 months vs 7.5 months). In this small cohort of patients we identified a higher than expected rates of KMT2C mutations. A KMT2C mutation was identified in 67% of samples which is higher than previously described rates of 7-15% in TCGA and METABRIC cohorts respectively. The presence of a concordant ctDNA KMT2C mutation in baseline and progression samples was a poor prognostic indicator (mPFS 7.5 months) when compared with no KMT2C mutation (mPFS 15.5 months) p =0.14 suggesting a possible mechanism of resistance to tucatinib. Conclusions: In this retrospective analysis of patients treated with Tucatinib, Trastuzumab and Capecitabine we established the genomic landscape of heavily pre-treated, advanced HER2-positive disease is complex. We observed there may be a correlation with KMT2C mutations and poorer clinical outcomes, however given the very small number of patients this data remains hypothesis generating, requiring validation with larger datasets. Citation Format: Elizabeth Blackley, Courtney van Geelen, Yi-An Ko, Stephen Wong, Miriam Yeung, Stephen Luen, Sarah-Jane Dawson, Sherene Loi. Investigation of the genomic evolution of HER2-positive breast cancer following progression on dual HER2-targetted therapy with Trastuzumab and Tucatinib [abstract]. In: Proceedings of the 2023 San Antonio Breast Cancer Symposium; 2023 Dec 5-9; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2024;84(9 Suppl):Abstract nr PO2-04-06.

Genome-wide Expanding of Genetic Evolution and Potential Pathogenicity in Vibrio alginolyticus

Vibrio alginolyticus, an emergent species of Vibrio genus, exists in aquatic and marine environments. It has undergone genetic diversification, but its detailed genomic diversity is still unclear. Here, we performed a multi-dimensional comparative genomic analysis to explore the population phylogeny, virulence-related genes and potential drug resistance genes of 184 V. alginolyticus isolates. Although genetic diversity is complex, we analyzed the population structure using three sub-datasets, including the subdivision for three lineages into sublineages and the distribution of strains in the marine ecological niche. Accessory genes, most of which reclassified V. alginolyticus genomes as different but with relatively close affinities, were nonuniformly distributed among these isolates. We demonstrated that the spread of some post-evolutionary isolates (mainly L3 strains isolated from Chinese territorial seas) was likely to be closely related to human activities, whereas other more ancestral strains (strains in the L1 and L2) tended to be locally endemic and formed clonal complex groups. In terms of pathogenicity, the potential virulence factors were mainly associated with toxin, adherence, motility, chemotaxis, and the type III secretion system (T3SS). We also found five types of antibacterial drug resistance genes. The prevalence of β-lactam resistance genes was 100%, which indicated that there may be a potential risk of natural resistance to β-lactam drugs. Our study reveals insights into genomic characteristics, evolution and potential virulence-associated gene profiles of V. alginolyticus.

Single-cell RNA-seq of Drosophila miranda testis reveals the evolution and trajectory of germline sex chromosome regulation.

Although sex chromosomes have evolved from autosomes, they often have unusual regulatory regimes that are sex- and cell-type-specific such as dosage compensation (DC) and meiotic sex chromosome inactivation (MSCI). The molecular mechanisms and evolutionary forces driving these unique transcriptional programs are critical for genome evolution but have been, in the case of MSCI in Drosophila, subject to continuous debate. Here, we take advantage of the younger sex chromosomes in D. miranda (XR and the neo-X) to infer how former autosomes acquire sex-chromosome-specific regulatory programs using single-cell and bulk RNA sequencing and ribosome profiling, in a comparative evolutionary context. We show that contrary to mammals and worms, the X down-regulation through germline progression is most consistent with the shutdown of DC instead of MSCI, resulting in half gene dosage at the end of meiosis for all 3 X's. Moreover, lowly expressed germline and meiotic genes on the neo-X are ancestrally lowly expressed, instead of acquired suppression after sex linkage. For the young neo-X, DC is incomplete across all tissue and cell types and this dosage imbalance is rescued by contributions from Y-linked gametologs which produce transcripts that are translated to compensate both gene and protein dosage. We find an excess of previously autosomal testis genes becoming Y-specific, showing that the neo-Y and its masculinization likely resolve sexual antagonism. Multicopy neo-sex genes are predominantly expressed during meiotic stages of spermatogenesis, consistent with their amplification being driven to interfere with mendelian segregation. Altogether, this study reveals germline regulation of evolving sex chromosomes and elucidates the consequences these unique regulatory mechanisms have on the evolution of sex chromosome architecture.

Insights into chloroplast genome evolution in Rutaceae through population genomics

Chloroplast genomes, pivotal for understanding plant evolution, remain unexplored in Rutaceae, a family with key perennial crops like citrus. Leveraging next-generation sequencing data from 509 Rutaceae accessions across 15 species, we conducted a de novo assembly of 343 chloroplast genomes, unveiling a chloroplast variation map highlighting the heterogeneous evolution rates across genome regions. Notably, differences in chloroplast genome size primarily originate from large single-copy and small single-copy regions. Structural variants predominantly occurred in the single-copy region, with two insertions located at the single-copy and inverted repeat region boundary. Phylogenetic analysis, principal component analysis, and population genetic statistics confirmed the cohesive clustering of different Citrus species, reflecting evolutionary dynamics in Citrus diversification. Furthermore, a close chloroplast genetic affinity was revealed among Atalantia (previously regarded as primitive citrus), Clausena, and Murraya. Zanthoxylum formed a distinct group with heightened genetic diversity. Through expanding our analysis to include 34 published chloroplast genomes, we explored chloroplast gene selection, revealing divergent evolutionary trends in photosynthetic pathways. While Photosystem I and Photosystem II exhibited robust negative selection, indicating stability, the Nicotinamide adenine dinucleotide (NADH) dehydrogenase pathway demonstrated rapid evolution, which was indicative of environmental adaptation. Finally, we discussed the effects of gene length and GC content on chloroplast gene evolution. In conclusion, our study reveals the genetic characterization of chloroplast genomes during Rutaceae diversification, providing insights into the evolutionary history of this family.

Growth and Genome Features of Non-O1/O139 Vibrio cholerae Isolated from Three Species of Common Freshwater Fish

Vibrio cholerae is the etiological agent of cholera in humans. The bacterium is frequently detected in aquatic products worldwide. However, the current literature on the genome evolution of V. cholerae of aquatic animal origins is limited. Here, we firstly characterized the growth and genome features of V. cholerae isolates with different resistance phenotypes from three species of common freshwater fish. The results revealed that the non-O1/O139 V. cholerae isolates (n = 4) were halophilic and grew optimally at 2% NaCl and pH 8.0. Their draft genome sequences were 3.89 Mb–4.15 Mb with an average GC content of 47.35–47.63%. Approximately 3366–3561 genes were predicted to encode proteins, but 14.9–17.3% of them were of an unknown function. A number of strain-specific genes (n = 221–311) were found in the four V. cholerae isolates, 3 of which belonged to none of any of the known sequence types (STs). Several putative mobile genetic elements (MGEs) existed in the V. cholerae isolates, including genomic islands (n = 4–9), prophages (n = 0–3), integrons (n = 1–1), and insertion sequences (n = 0–3). Notably, CRISPR-Cas system arrays (n = 2–10) were found in the V. cholerae genomes, whereby the potential immunity defense system could be active. Comparative genomic analyses also revealed many putative virulence-associated genes (n = 106–122) and antibiotic resistance-related genes (n = 6–9). Overall, the results of this study demonstrate the bacterial broader-spectrum growth traits and fill prior gaps in the genomes of V. cholerae originating from freshwater fish.

Discovery and genome-wide characterization of a novel miniature inverted repeat transposable element reveal genome-specific distribution in Glycine.

Miniature inverted repeat transposable elements (MITEs) are a dynamic component responsible for genome evolution. Tourist MITEs are derived from and mobilized by elements from the harbinger superfamily. In this study, a novel family of Tourist-like MITE was characterized in wild soybean species Glycine falcata. The new GftoMITE1 was initially discovered as an insertional polymorphism of the centromere-specific histone H3 (CenH3) gene in G. falcata. Using polymerase chain reaction, cloning and sequencing approaches, we showed a high number of copies of the GftoMITE1 family. Extensive bioinformatic analyses revealed the genome-level distribution and locus-specific mapping of GftoMITE1 members in Glycine species. Our results provide the first extensive characterization of the GftoMITE1 family and contribute to the understanding of the evolution of MITEs in the Glycine genus. Genome-specific GftoMITE1 was prominent in perennial wild soybean species, but not in annual cultivated soybean (Glycine max) or its progenitor (Glycine soja). We discuss that the GftoMITE1 family reveals a single rapid amplification in G. falcata and could have potential implications for gene regulation and soybean breeding as an efficient genetic marker for germplasm utilization in the future.

Transposable elements in Rosaceae: insights into genome evolution, expression dynamics, and Syntenic gene regulation

Abstract Transposable elements (TEs) exert significant influence on plant genomic structure and gene expression. Here, we explored TE-related aspects across fourteen Rosaceae genomes, investigating genomic distribution, transposition activity, expression patterns, and nearby differentially expressed genes (DEGs). Analyses unveiled distinct long terminal repeat retrotransposon (LTR-RT) evolutionary patterns, reflecting varied genome size changes among nine species over the past million years. In the past 2.5 million years, Rubus idaeus showed a transposition rate twice as fast as Fragaria vesca, while Pyrus bretschneideri displayed significantly faster transposition compared to Crataegus pinnatifida. Genes adjacent to recent TE insertions were linked to adversity resistance, while those near previous insertions were functionally enriched in morphogenesis, enzyme activity, and metabolic processes. Expression analysis revealed diverse responses of LTR-RTs to internal or external conditions. Furthermore, we identified 3695 pairs of syntenic DEGs proximal to TEs in Malus domestica cv. Gala and Malus domestica (GDDH13), suggesting TE insertions may contribute to varietal trait differences in these apple varieties. Our study across representative Rosaceae species underscores the pivotal role of TEs in plant genome evolution within this diverse family. It elucidates how these elements regulate syntenic DEGs on a genome-wide scale, offering insights into Rosaceae-specific genomic evolution.

The genome of the Arctic snow-alga Limnomonas spitsbergensis (Chlamydomonadales).

Snow-algae are a diverse group of extremophilic microeukaryotes found on melting polar and alpine snowfields. They play an important role in the microbial ecology of the cryosphere, and their propagation on snow and ice surfaces may in part accelerate climate-induced melting of these systems. High quality snow-algae genomes are needed for studies on their unique physiology, adaptive mechanisms and genome evolution under multiple forms of stress, including cold temperatures and intense sunlight. Here we assembled and annotated the genome of Limnomonas spitsbergensis, a cryophilic biciliate green alga originally isolated from melting snow on Svalbard, in the Arctic. The L. spitsbergensis genome assembly is based primarily on the use of PacBio long reads and secondly Illumina short reads, with an assembly size of 260.248 Mb in 124 contigs. A combination of three alternative annotation strategies were used including protein homology, RNA-seq evidence and PacBio full length transcript isoforms. The best merged set of annotations identified 18,277 protein-coding genes, which were 95.2% complete based on BUSCO analysis. We also provide the annotated mitogenome, which is a relatively large 77.942 kb circular mapping sequence containing extensive repeats. The L. spitsbergensis genome will provide a new resource for research on snow-algae adaptation, behavior and natural selection in unique, low-temperature terrestrial environments that are under threat from climate change.

Shared features underlying compact genomes and extreme habitat use in chironomid midges.

Non-biting midges (family Chironomidae) are found throughout the world in a diverse array of aquatic and terrestrial habitats, can often tolerate harsh conditions such as hypoxia or desiccation, and have consistently compact genomes. Yet we know little about the shared molecular basis for these attributes and how they have evolved across the family. Here, we address these questions by first creating high-quality, annotated reference assemblies for Tanytarsus gracilentus (subfamily Chironominae, tribe Tanytarsini) and Parochlus steinenii (subfamily Podonominae). Using these and other publicly available assemblies, we created a time-calibrated phylogenomic tree for family Chironomidae with outgroups from order Diptera. We used this phylogeny to test for features associated with compact genomes, as well as examining patterns of gene family evolution and positive selection that may underlie chironomid habitat tolerances. Our results suggest that compact genomes evolved in the common ancestor of Chironomidae and Ceratopogonidae, and that this occurred mainly through reductions in non-coding regions (introns, intergenic sequences, and repeat elements). Significantly expanded gene families in Chironomidae included biological processes that may relate to tolerance of stressful environments, such as temperature homeostasis, carbohydrate transport, melanization defense response, and trehalose transport. We identified several positively selected genes in Chironomidae, notably sulfonylurea receptor, CREB-binding protein, and protein kinase D. Our results improve our understanding of the evolution of small genomes and extreme habitat use in this widely distributed group.

Surfing the Waves of SARS-CoV-2: Analysis of Viral Genome Variants Using an NGS Survey in Verona, Italy

The availability of new technologies for deep sequencing, including next-generation sequencing (NGS), allows for the detection of viral genome variations. The epidemiological determination of SARS-CoV-2 viral genome changes during the pandemic waves displayed the genome evolution and subsequent onset of variants over time. These variants were often associated with a different impact on viral transmission and disease severity. We investigated, in a retrospective study, the trend of SARS-CoV-2-positive samples collected from the start of the Italian pandemic (January 2020) to June 2023. In addition, viral RNAs extracted from 938 nasopharyngeal swab samples were analyzed using NGS between February 2022 and June 2023. Sequences were analyzed with bioinformatic tools to identify lineages and mutations and for phylogenetic studies. Six pandemic waves were detected. In our samples, we predominantly detected BA.2, BQ.1, BA.5.1, BA.5.2, and, more recently, XBB.1 and its subvariants. The data describe the SARS-CoV-2 genome evolution involved in viral interactions with the host and the dynamics of specific genome mutations and deletions.

Long-read sequencing of extrachromosomal circular DNA and genome assembly of a Solanum lycopersicum breeding line revealed active LTR retrotransposons originating from S. Peruvianum L. introgressions

Transposable elements (TEs) are a major force in the evolution of plant genomes. Differences in the transposition activities and landscapes of TEs can vary substantially, even in closely related species. Interspecific hybridization, a widely employed technique in tomato breeding, results in the creation of novel combinations of TEs from distinct species. The implications of this process for TE transposition activity have not been studied in modern cultivars. In this study, we used nanopore sequencing of extrachromosomal circular DNA (eccDNA) and identified two highly active Ty1/Copia LTR retrotransposon families of tomato (Solanum lycopersicum), called Salsa and Ketchup. Elements of these families produce thousands of eccDNAs under controlled conditions and epigenetic stress. EccDNA sequence analysis revealed that the major parts of eccDNA produced by Ketchup and Salsa exhibited low similarity to the S. lycopersicum genomic sequence. To trace the origin of these TEs, whole-genome nanopore sequencing and de novo genome assembly were performed. We found that these TEs occurred in a tomato breeding line via interspecific introgression from S. peruvianum. Our findings collectively show that interspecific introgressions can contribute to both genetic and phenotypic diversity not only by introducing novel genetic variants, but also by importing active transposable elements from other species.

Emx2 underlies the development and evolution of marsupial gliding membranes

Phenotypic variation among species is a product of evolutionary changes to developmental programs1,2. However, how these changes generate novel morphological traits remains largely unclear. Here we studied the genomic and developmental basis of the mammalian gliding membrane, or patagium—an adaptative trait that has repeatedly evolved in different lineages, including in closely related marsupial species. Through comparative genomic analysis of 15 marsupial genomes, both from gliding and non-gliding species, we find that the Emx2 locus experienced lineage-specific patterns of accelerated cis-regulatory evolution in gliding species. By combining epigenomics, transcriptomics and in-pouch marsupial transgenics, we show that Emx2 is a critical upstream regulator of patagium development. Moreover, we identify different cis-regulatory elements that may be responsible for driving increased Emx2 expression levels in gliding species. Lastly, using mouse functional experiments, we find evidence that Emx2 expression patterns in gliders may have been modified from a pre-existing program found in all mammals. Together, our results suggest that patagia repeatedly originated through a process of convergent genomic evolution, whereby regulation of Emx2 was altered by distinct cis-regulatory elements in independently evolved species. Thus, different regulatory elements targeting the same key developmental gene may constitute an effective strategy by which natural selection has harnessed regulatory evolution in marsupial genomes to generate phenotypic novelty.

Massive genome reduction predates the divergence of Symbiodiniaceae dinoflagellates.

Dinoflagellates in the family Symbiodiniaceae are taxonomically diverse, predominantly symbiotic lineages that are well-known for their association with corals. The ancestor of these taxa is believed to have been free-living. The establishment of symbiosis (i.e., symbiogenesis) is hypothesised to have occurred multiple times during Symbiodiniaceae evolution, but its impact on genome evolution of these taxa is largely unknown. Among Symbiodiniaceae, the genus Effrenium is a free-living lineage that is phylogenetically positioned between two robustly supported groups of genera within which symbiotic taxa have emerged. The apparent lack of symbiogenesis in Effrenium suggests that the ancestral features of Symbiodiniaceae may have been retained in this lineage. Here we present de novo assembled genomes (1.2-1.9 Gbp in size) and transcriptome data from three isolates of Effrenium voratum and conduct a comparative analysis that includes 16 Symbiodiniaceae taxa and the other dinoflagellates. Surprisingly, we find that genome reduction, which is often associated with a symbiotic lifestyle, predates the origin of Symbiodiniaceae. The free-living lifestyle distinguishes Effrenium from symbiotic Symbiodiniaceae vis-à-vis their longer introns, more-extensive mRNA editing, fewer (~30%) lineage-specific gene sets, and lower (~10%) level of pseudogenisation. These results demonstrate how genome reduction and the adaptation to distinct lifestyles intersect to drive diversification and genome evolution of Symbiodiniaceae.